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Owing to the variety of possible charge and spin states and to the different ways of coupUng to the 
environment, paramagnetic centres in wide band-gap semiconductors and insulators exhibit a strikingly 
rich spectrum of properties and functionalities, exploited in commercial light emitters and proposed for 
applications in quantum information. Here we demonstrate, by combining synchrotron techniques with 
magnetic, optical and ab initio studies, that the codoping of GaNiMn with Mg allows to control the Mn"^ 
charge and spin state in the range 3<i3<5 and 2>5>1. According to our results, this outstanding degree of 
tunability arises from the formation of hitherto concealed cation complexes Mn-Mg/^ where the number of 
ligands icis pre-defined by fabrication conditions. The properties of these complexes allow to extend towards 
the infrared the already remarkable optical capabilities of nitrides, open to solotronics functionalities, and 
generally represent a fresh perspective for magnetic semiconductors. 



Group III nitrides currently dominate the field of visible and ultraviolet photonics, due to the flexibility 
of gap engineering, to the availability of both n- and p-type material, to their high thermal stability, 
and large heat conductivity (ref. 1 and references therein). Moreover, recent studies of magnetically 
doped semiconductors^"^ and oxides^'^ have pointed out how crucial is the spatial distribution of dopants in 
the host matrix. In particular, the aggregation of magnetic impurities resulting at the nanoscale in the 
formation of chemical and/or crystallographic inhomogeneities has a dramatic influence on the magnetic, 
electronic and optical properties of the systems. The possibility to reach a control over the formation, 
structure, arrangement, and effect of inhomogeneities on the chemical and physical behavior of the host 
material, can be envisaged to have outstanding consequences in the design and functionality of the next 
generation of devices. The manipulation of the charge-state of paramagnetic centres in wide band gap 
semiconductors by modifj^ing the local electric field with a gate voltage^ or chemically by changing the 
termination of the surface^, has already been reported. 

Here, we demonstrate the efficiency of an alternative approach based on the controlled formation of com- 
plexes, involving one magnetic impurity and one or more electrically active dopant. Peculiar arrangements in the 
electronic structure of simple paramagnetic complexes have previously been theoretically investigated^". With 
synchrotron radiation methods supported by ab initio computations, we show that codoping of GaN with Mn and 
Mg results in the formation of cation complexes Mn-Mg;^. Depending on the number of ligands k predefined by 
fabrication conditions, a strong electron-phonon interaction and consequently an efficient and broadband 
infrared (IR) photoluminescence (PL) are promoted. Furthermore, we find that the Mn-Mg^ cation complexes 
allow to control the charge and spin states of the transition metal (TM) ions, with prospective implications in the 
design of the coupling between localized spins in magnetic semiconductors and in the optimisation of centres for 
solotronic applications ^\ 
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Results 

The samples studied are single- crystal GaN layers, codoped with Mn 
and Mg, each with concentrations lower than 1%. Their fabrication, 
architecture, and preliminary characterisation are summarised in the 
Methods. Aware that the incorporation of interstitial hydrogen form- 
ing complexes with Mg is a recurrent challenge associated with Mg- 
doping of GaN, we have performed a careful chemical analysis, 
detailed in the Supplementary Figure SI, showing that in our case 
we can reasonably neglect the effect of interstitial hydrogen. 

Nature of Mn-Mg complexes. Having established the single 
crystallinity and chemical homogeneity of the samples via a 
spectrum of both local and averaging characterisation techniques, 
we discuss the lattice positions of Mn and Mg impurities in GaN, as 
determined from extended x-ray absorption fine structure (EXAFS) 
measurements at the Mn K-edge and by ab initio computations. We 
find that within a confidence of 90% the Mn ions occupy exclusively 
Ga-substitutional positions. This result is consistent with the ah 
initio computations reported in the Supplementary Table SI 
showing that the energy required to place a Mn or Mg ion in 
either octahedral or tetrahedral interstitial positions of GaN is 
more than 4 eV higher than the one needed to incorporate it in a 
Ga-substitutional site. 

According to our previous EXAFS and electron energy loss spec- 
troscopy (EELS) results, Mn is randomly and homogeneously dis- 
tributed in GaN:Mn, at least up to a Mn concentration of 3% (refs. 
12,13). While in the conventional treatment of dilute magnetic semi- 
conductors, the spatial distributions of co-dopants and TM ions are 
assumed to be uncorrected, a quantitative analysis of our EXAFS 
data points to a substantial correlation between the positions occu- 
pied by Mn and Mg in the host lattice. Simulations of EXAFS spectra 
for a large variety of relaxed defects in GaN were performed and 
reported in the Supplementary Information, and indicate that the 
combination of substitutional Mn and Mg is the most likely to 
account for the experimental EXAFS data. The experimental 
EXAFS spectra are then fitted, according to the procedure described 
in the Supplementary Information. As shown in Fig. lb, the number 
of substitutional Mg atoms (/iMg) in the first cation coordination 



sphere of Mn increases linearly with the ratio between the Mg and 
Mn concentrations, y = Xug/xun^ up to y = 3 and then saturates at 
higher;/ values. Simultaneously, the average distance beween Mn and 
the nearest neighbour N atoms dun-N diminishes in range up to y = 
3 and then levels off, as seen in Fig. le. The experimental EXAFS 
spectra from which (^Mn-N and are extracted are plotted in 
Supplementary Figure S4 for different values of y. It is possible to 
confirm the correlation between the Mn and Mg positions by com- 
paring the spectra for y > 0 (correlated) and y = 0, the latter being 
strictly equivalent to the not correlated case where Mn does not 
interact with any Mg atom. 

These results can be explained by our ah initio computations. In 
particular, the estimated values of the pairing energies shown in 
Fig. la for up to three Mg per Mn are all negative, demonstrating the 
tendency of Mg to form the Mn-Mg;^ complexes sketched in Fig. la 
and in Supplementary Figure S8. Taking into account the statistical 
distribution of k, implying that some of these complexes can coexist 
at a given y, we have obtained with no adjustable parameters a 
remarkable agreement between the experimental and computed 
trends, describing how the number of bound Mg atoms (Fig. lb) 
and the shortening of the bond length (Fig. lc,d,e) depend on y. 
This agreement with the theory based on the statistical distribution 
of the complexes populations detailed in the Methods implies a com- 
parable stability of the complexes with different k, wp to k = 3. 
According to the DFT calculations and as reported in Fig. la, a 
similar binding energy (0.5 eV) is indeed computed for Mn-Mg 
and Mn-Mg2 but a lower value (0.2 eV) is expected for Mn-Mg3. If 
these DFT predictions are quantitatively correct, we may await a 
concentration of Mn-Mg3 complexes somewhat lower than the one 
resulting from Psiy) for y < 1. According to the ah initio results in 
Fig. Id, the bond shortening is particularly significant for the Mn-N 
pairs nearest to Mg. Intuitively, this effect results from a sizable 
charge transfer at the Mn-Mg;^ complexes, that leads to the onset 
of Coulomb attraction between the ionized Mg acceptors and Mn 
donor. 

Control of charge and spin state. Having discussed the atomistic 
structure of the Mn-Mg^ complexes, we determine their charge and 
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Figure 1 | Mn-Mg^complexes predicted by theory and experimental demonstration by EXAFS. (a) most stable Mn-Mg^t complexes (/c= 1,2 and 3) and 
their pairing energies - computed by DFT within the GGA+U approximation - relatively to the previous Mn-Mgjt_i complex. A gray arrow indicates the 
[0001] direction, i.e. the GaN c-axis; (b) number of Mg atoms seen by Mn in the first cation coordination sphere (/^Mg)) extracted from EXAFS 
measurements and DFT predictions; (c) schematic representation of a Mn-Mg complex; (d) detailed values of the bond lengths for different complexes, 
calculated by DFT; (e) average Mn-N bond length (^^mii-n) as a function of the Mg/Mn ratio, from EXAFS analysis and DFT calculations. 
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spin states. According to our ab initio computations, if Mn-Mg/^ 
complexes were not formed, only a single hole would be trapped 
by Mn upon Mg codoping, resulting in the Mn spin state S =3/2. 
In that case, \i y > 1, extra holes would be directed to the valence 
band, increasing thep-type conductivity. However, layers containing 
a small concentration of Mg (y < 3) are insulating and the onset of 
conductivity takes place only iox y > 3, with a decrease in the 
resistance of the layers of more than three orders of magnitude, 
indicating that at lowy, Mg are bound to Mn in complexes, and do 
not exist as isolated acceptors, as free Mg acceptors are only present 
for )/ > 3. Considering the Mn-Mg^ complexes, our computations 
predict the presence of Mn with spin state down to S = 1 for A; > 2. 
Altogether, our theoretical studies suggest the possibility to control, 
with Mg codoping, the Mn"^ charge and spin state over the range 3 < 
n < 5 and 2 > S > 1, respectively. 

We give now experimental verification of these theoretical expec- 
tations by tracing the evolution of S upon Mg codoping. We exploit 
the Kj5 x-ray emission spectroscopy (XES) technique which, by prob- 
ing the 3p ^ Is transitions, is sensitive to the magnitude of the 
exchange interaction between the Mn 3p core-hole and the net mag- 
netic moment in the Mn 3d valence shelP^. We note that compared to 
x-ray absorption spectroscopy (XAS), where a charge state depend- 
ent shift of the absorption edge is also visible, the XES data are more 
linearly correlated with the spin state, and depend less on the atomic 
configuration^^. The spectra are plotted in Fig. 2a from which we 
extract, by comparison with reference oxides, as described in the 
Supplementary Information, the values of S reported in Fig. 2b. 
This analysis demonstrates that upon Mg codoping the Mn spin state 
first decreases linearly from 2 to 1, and then saturates for ;/ > 3. As 
shown in Fig. 2b, these results are in accord with the magnitudes of S 
determined from the ah initio computations for a unit cell containing 
one Mn. 

Magnetism of complexes. Having established the TM charge and 
spin states, one can examine the properties of the centres and their 



coupling to the environment in terms of the time-honoured crystal- 
field theory^^, providing the expected structure and symmetries of the 
relevant energy levels at a given oxidation state^^. In order to gain 
more understanding at the microscopic leve^^ it is instructive to 
consider the actual spin distribution around the Mn ions and its 
evolution with Mg codoping. As shown in Fig. 2b and in the 
Supplementary Table S2, our ah initio computations confirm that 
in the case of GaN:Mn the total magnetic moment m = 4 is built of 
nic ~ 4.5 on the d shell of the Mn ion and ~ —0.5 ^b on thep 
orbitals of the N ligands. This means, in agreement with the notion of 
strong coupling limit for GaN:Mn (ref 19), that the wave function of 
the hole provided by the Mn is equally spread between Mn and the 
neighbouring N ions. As reported in Fig. 2d, the delocalization of 
holes is enhanced ioxk= 1, suggesting that Mn-Mg complexes can 
mediate spin-dependent interactions between Mn spins. On the 
other hand, as evidenced in Fig. 2e,f, for A; = 2 and 3 the delocali- 
zation is quenched, and this result is important for understanding the 
optical data discussed later. 

The changes of the Mn spin state and Mn-N bonding with Mg 
codoping, evaluated by EXAFS, XES, and ah initio studies affect the 
magnetic properties of the system. Quite generally, the values of the 
magnetization M{H) and its anisotropy are determined by the rel- 
evant spin Hamiltonian, whose form, appropriate for a given S, is 
determined by crystal symmetry, including local strains, whereas the 
magnitudes of the spin Hamiltonian parameters provide infor- 
mation on the coupling to the environment (including p - d hybrid- 
ization) and on the strength of the spin-orbit interaction. In 
GaN:Mn, a sizable difference in the magnetization curves M{H) is 
observed for two orientations of the magnetic field H with respect to 
the wurtzite c-axis, H ± c and H \\ c (refs. 12,20). A quantitative 
analysis of such data allowed to obtain the values of the spin 
Hamiltonian parameters for this case, with S = 2 and wurtzite sym- 
j^g^jy 12,20 the case of GaN:(Mn,Mg), as seen in Fig. Ic-d, the 
distortions of the nitrogen tetrahedron and, in particular, the short- 
ening of the Mn-N bond induced by the presence of Mg, is much 
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Figure 2 | Evolution of the Mn spin state, (a) non-resonant XES for GaN:Mn samples with and without Mg; (b) evolution of the nominal spin value 
with the ratio between Mg and Mn concentration, extracted from the analysis of XES spectra and calculated via DFT; (c) computed spin polarization 
density [Ap = p(|) — p{[)] for Muca in wurtzite (wz)-GaN. The positive and negative spin polarizations are represented by violet and grey colours, 
respectively; (d-f) difference — pjt-i> between the spin polarizations of the Mn-Mgjt complexes, with /: = 1, 2 and 3. The red and blue colours represent 
positive and negative values, respectively. The blue colour in (d) indicates the enhanced delocalization of spin polarization in Mn-Mg, whereas the red 
colour in (e) and (f) points to a gradual shift of the spin density to Mn in Mn-Mg2 and then in Mn-Mg3. In all plots the contour value is set to 
0.005 corresponding to about 2 X 10~^ /^bA"^. 
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Figure 3 | Magnetism of Mn-Mgj^ complexes. Magnetic anisotropy 
energy density as a function of y measured by SQUID magnetometry, and 
calculated assuming that it originates from Mn without Mg in the first 
cation coordination sphere. Inset: normalized magnetization curves M-H 
of GaNiMn containing 0.4% of Mn without and with Mg codoping (left 
and right panel, respectively) measured at 1.85 K. 



more pronounced than the elongation of the Mn-N bond parallel to 
the c-axis. Accordingly, the local anisotropy will now be determined 
by the position of one or more Mg in the first cation coordination 
sphere. Since there is no preferential orientation for the occupation of 
any of the 12 nearest neighbour positions, the local anisotropy and, 
thus, the orientation of the easy axis will be randomly distributed 
over the Mn-Mg;^ complexes. Hence, the presence of Mg leads to the 
disappearance of anisotropy with respect to the c-axis. Indeed, as 
shown in Fig. 3, we observe that the magnitude A of the magnetic 
anisotropy decays to zero with increasing y. In agreement with the 
model, A{y) follows the probability Po{y) that no Mg is bound by Mn 
at a given y, A{y) = A(0)(1 — y/3)l At the same time, the absolute 
values of M{H) are determined by unknown parameters of the spin 
Hamiltonians, appropriate for particular values of S, local strains, 
and strength of the spin-orbit interaction. 

Effect of complexes on the luminescence of Mn. Particularly 
remarkable is the influence of the Mn-Mg^ complexes on the light 
emission of GaN. At present, TM-doping of oxides and II-VI 
chalcogenides allows the fabrication of optically pumped tunable 

lasers^^'^^. 

The PL spectra reported in Fig. 4a for our samples, extends for 
high y- ratios over a broad infrared band, covering two of the tele- 
communication windows, namely 1.33 ^m and 1.55 ^m. In addition, 
the multisite structure suggests that only one transition is involved, 
but different environments as well as phonon replicas are responsible 
for the peaks multiplicity. The overall shape is indicative of a Huang- 
Rhys factor Shr > 2. According to Fig. 4b, the PL intensity increases 
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Figure 4 | Infrared photoluminescence of GaN:(Mn,Mg) samples, (a) evolution of the PL spectra (excited with a 442 nm (2.8 eV) laser) as a function of 
the / ratio, measured at 2 K (a multiplying factor of 1.4 has been applied between consecutive spectra for clarity); (b) integrated PL intensity normalized 
by the Mn concentration and sample thickness as a function of the y ratio (red circles, left scale) and evolution of the fraction of different complexes, 
calculated as a function of / (lines, right scale); (c) evolution of the PL spectra with temperature, for a sample with y = 4. 1; the background at 296 K is also 
included. A factor of 1.4 has been applied between consecutive spectra for clarity, (d) evolution of the integrated PL intensity as a function of the inverse 
temperature, normalized to the value at T = 5 K; the line is a guide for the eye. Inset: levels relevant to PL as discussed in the Supplementary Information. 
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with the concentration of the Mn-Mg3 complexes Psiy) = (y/^Y, 
demonstrating that Mn^^, in presumably two slightly different envir- 
onments, accounts for the photoluminescence in GaN:(Mn,Mg). 

Importantly, as one sees in Fig 4c and d, the spectrally broad 
infrared emission persists up to room temperature, and is very 
attractive for ultrashort pulse generation as well as for wide infrared 
tunability. 

In order to understand the photoluminescence related to Mn-Mg;^ 
cation complexes, it is necessary to consider first the reason of the 
poor luminescence of Mn^^ in GaN. In contrast to Cr^^ in AI2O3 and 
ZnSe, Mn^^ in GaN does not show the strong and application-rel- 
evant red or infrared emission associated with optical transitions 
between the and crystal field multiplets of the 3d^ shell. This 
surprising result is however in agreement with previous extensive 
studies of GaN:Mn (refs 23-26), in which the intra-centre photolu- 
minescence of Mn^^ was found to be hardly detectable. This puzzling 
result can be traced back to an abnormally small value of the Huang- 
Rhys factor Shr < 1 implied by optical absorption^^'^^, this obser- 
vation being interpreted^^ in terms of a strong p-d coupling in 
GaN:Mn, leading to a significant admixture of anion t2 orbitals with 
the wave functions of the and ^T2 states. Such delocalization of the 
centre wave function, clearly visible in Fig. 2c, reduces the electron- 
phonon coupling. As a result, the oscillator strength is shifted to the 
zero-phonon line at the expense of the phonon-assisted transitions 
that cease to constitute the channel for PL excitation. 

On the other hand, highly efficient and Stokes shifted PL, similar 
to the one presented in Fig. 4, was found for GaN:(Mn,Mg) and 
assigned to intra Mn^^ transitions^^'^^'^^. However, we have demon- 
strated here that following the evolution of the PL as a function of the 
different populations of complexes (Fig. 4b), the Mn^^ present in the 
complexes are likely to be at the origin of the infrared PL signal. This 
assignment is consistent with the shape of the PL spectrum in a 10 T 
magnetic field - given in the Supplementary Figure S7 - which can be 
described by a splitting of the ground state into three components. 
This splitting is expected for S = 1, as observed previously^^ for V^^ 
in AIN. 

In order to explain the origin of the PL activation upon codoping 
with Mg evidenced in Fig. 4a, we refer to Fig. 2d-f where the varia- 
tions in the local spin densities brought about by the binding of an 
increasing number of Mg ions are shown. According to these data, 
the complexing with one Mg ion enhances the delocalization of the 
spin density over neighbouring N anions. However, with the binding 
of two and then three Mg ions, the delocalization of the spin density 
decreases. Accordingly, a strong electron-phonon coupling and, 
thus, a large magnitude of Shr is restored, particularly for Mn- 
Mg3 complexes, where two d electrons reside in the e orbitals of 
the ^A2 ground state. 

This IR broadband PL promoted by Mn-Mg/^ cation complexes is 
of high relevance in e.g. laser and telecommunication technologies. 
Actually, in comparison with AI2O3 and ZnSe, GaN has respectively 
seven and twelve times better thermal conductivity, lessening ther- 
mal effects even for high laser powers and intensities. 

Discussion 

We conclude that the data presented here indicate a new way to 
manipulate the charge and spin state of single paramagnetic centres 
by complexing a magnetic impurity with electrical dopants. The 
demonstration of these new degrees of freedom opens wide prospects 
illustrated here by the infrared emission of GaN:(Mn,Mg). Another 
line of research is to explore the potential of these Mn-Mg^ cation 
complexes for mediating the coupling between localized spins in 
magnetic semiconductors. Furthermore, these centres may serve 
for storing and manipulating information in a single qubit or for 
single photon generation ^\ Interestingly, unlike the case of 
CdTe:Mn (ref 31) or InAs:Mn (ref 32), it is not necessary to place 
GaN:(Mn,Mg) in a quantum dot, as Mn in GaN can bind the 



exciton^^ needed to read or write information. The possibility to 
change energy level splitting, excitation channel and excited state 
lifetime by manipulating the Mn charge and spin state through 
Mn-Mg^ cation complexes offers a not yet explored spectrum of 
opportunities for further investigations. 

Methods 

Growth and preliminary characterisation. The samples consist of single crystal 
wurtzite (wz) GaN codoped with Mn and Mg grown by metalorganic vapor phase 
epitaxy (MOVPE) on a 1 fim GaN buffer layer on c-plane sapphire, according to the 
procedure described elsewhere^^'^^ The doped layer is 600 nm thick. The samples are 
grown under H2 atmosphere, with a pressure of 200 mbar and a temperature of 
850°C. The precursors used are ammonia (NH3) for nitrogen, trimethylgallium 
(TMGa) for Ga, dicyclopentadienyl-magnesium (Cp2Mg) for Mg and 
dicyclopentadienyl- manganese (Cp2Mn) for Mn. The source flow of ammonia was 
kept constant at 1500 seem, the TMGa at 5 seem, and Cp2Mg was varied between 150 
and 450 seem as Cp2Mn was varied between 75 and 490 seem. The Mn and Mg 
concentrations considered in this work are both between 0 and 1% as measured by 
secondary ion mass spectroscopy (SIMS). The absence of parasitic elements like 
hydrogen or oxygen has been carefully checked with SIMS, energy dispersive x-ray 
spectropscopy (EDX), Raman spectropscopy, and electron energy loss spectroscopy 
(EELS). Prior to the extensive synchrotron investigations by EXAFS and XES, the 
structure of the layers has been characterised by high -resolution x-ray diffraction 
(HRXRD) on a X'Pert PRO MRD system with a dynamics as high as 10^ between the 
GaN (002) peak and the noise. In addition, high-resolution transmission electron 
microscopy (HRTEM) was performed on a JEOL 2011 Fast TEM microscope 
operating at 200 kV and capable of an ultimate point-to-point resolution of 0.19 nm 
and allowing to image lattice fringes with a 0.14-nm resolution. The combination of 
the two techniques has allowed us to rule out the presence of precipitation in the 
layers. 

EXAFS. EXAFS spectroscopy has been carried out at the BM08-GILDA Italian 
beamline^"^ at the ESRF (Grenoble, France). The Mn K edge x-ray absorption spectra 
have been acquired using a monochromator equipped with a pair of Si(311) crystals 
and run in dynamical focusing mode. Harmonics rejection is achieved through a pair 
of Pd-coated mirrors and the monochromator de-tuning. The data are collected in the 
fluorescence mode using a 13 -element hyperpure Ge detector and normalized by the 
incoming flux measured with an ion chamber. The incident beam is at 55.7° in respect 
to the sample surface to avoid dichroic effects. The samples are cooled down to liquid 
nitrogen temperature. The counting time and the number of scans for each sample 
have been chosen in order to collect at least 10^ counts per point. The EXAFS signal, 
Xik), is extracted from the absorption data, fi{E), using a smoothing spline algorithm 
(as implemented in the viper program) and choosing the energy edge, Eq, at the 
maximum of the derivative. The data analysis is detailed in the Supplementary 
Information. 

XES. Non resonant XES has been measured at the ID26 beamline of the ESRF^^. The 
optics for the incoming beam consists of three coupled undulators, a double Si crystal 
monochromator and three Si coated mirrors working at 2.5 mrad incidence for 
harmonics rejection and beam focusing. The emission spectrometer is run in a 
vertical Rowland geometry with five Si(llO) analyzer crystals working at the (440) 
reflection, that is, around a Bragg angle of 84.2°. The spectrometer energy broadening 
is approximately 0.9 eV at Mn K-edge (6539 eV). The experimental geometry 
consists of the spectrometer and the incoming beam at 90° on the same scattering 
plane (to minimize the elastic contribution) with the sample surface placed vertically 
to this plane and at 55.7° incidence angle, that is, the magic angle for wurtzite 
symmetry in order to avoid dichroism effects due to the linearly polarized beam^*^. All 
the samples are measured at room temperature and are tested against radiation 
damage. The emitted fluorescence is scanned around the Mn KjS main line with the 
incoming excitation at 6700 eV. The quantitative data analysis is based on the 
integrated absolute values of the difference spectra (IAD) and is performed as a 
function of the Mg/Mn concentration ratio, y. The data have been calibrated by the 
IAD values obtained from commercial Mn-oxides powders, assuming the ionic 
approximation and considering the high-spin scenario for both systems. The details 
of the analysis are reported in the Supplementary Information. 

SQUID magnetometry. The magnetic anisotropy energy density has been assessed 
by integrating the area between the magnetization curves measured along easy {H ± 
c) and hard (H 1 1 c) directions and expressed as energy per one Mn atom in the given 
layer. Magnetization curves are measured at 1.85 K using a superconductor quantum 
interference device (SQUID), as described previously^^'^^. The size of the error bars for 
the anisotropy is determined mostly by the errors related to the inaccuracy of 
substrate signal compensation and do not include the uncertainty generated by the 
insufficient strength of the magnetic field to saturate M for hard direction in our 
SQUID magnetometer (50 kOe). 

PL. Photoluminescence is excited with a continuous wave 404 nm (3.1 eV) or 
442 nm (2.8 eV) laser with the excitation power up to tens of mW. The temperature 
has been varied in the range between 2 K and 296 K. An InGaAs type CCD camera 
coupled to a grating (either 300 grooves/mm or 1200 grooves/mm) monochromator 
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is used as detector. A long wavelength pass filter is placed at the entrance of the 
monochromator for cutting off the stray laser light. The detection is carried out in the 
range from 0.7 eV to 1.5 eV with a spectral resolution of 0.5 meV. The integration of 
the PL signal in Fig. 4b and 4d has been performed between 900 meV and 1 100 meV. 
The magnetooptical measurements reported in the Supplementary Information are 
performed in Faraday configuration (5 1 1 A;) using a cryostat equipped with a 
superconducting coil providing a magnetic field up to 10 T. 

Theory - DFT. Calculations for Mn-Mg^ complexes in wz-GaN are performed 
within the GGA+U approximation using the Quantum Espresso code^^. A 96-atoms 
supercell and a 3X3X3 Monkhorst-Pack grid for Brillouin zone sampling are 
employed. The pairing energies are calculated from the following formula: AE = 
£tot(MnGa47_fcN48 : Mg^ + Ga48N48) - £tot(MnGa47_^+iN48 : Mgfc_i + Ga47N48 : Mg), 
for k ranging from 1 to 5. From the computed values of the magnetic moment and 
from the Mn-N distance of every possible complexes, one can obtain, taking into 
account the relative statistical weight Pkiy) of particular Mn-Mg^ configurations, the 
average variation of these values as a function of the Mg/Mn ratio. The values of these 
Pkiy) are then approximated with a binomial law, considering that if Mn can bind up 
to m Mg atoms (m = 3 in our case), the occurrence probability Pkiy) of particular 
complexes Mn-Mg^, at a given ratio of the Mg to Mn concentration 7, is given by the 
binomial distribution, Pkiy) = {'k)iy/^)^i^ —y/m)^~^ for)/ < m, whereas for)/ > 
m, Pmiy) = 1 and Pkiy) = 0 for k < m. The computed values of the local magnetic 
moment on Mnoa and its nearest neighbouring N atoms as well as its magnitude in 
the Mn unit cell are collected in the Supplementary Table S2, whereas the contour plot 
of the spin polarization is shown in Fig. 2c. 
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